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TITLE: 

TYPES 

NUMBER OF WORDS: 

TEMPORARY STORAGE: 

DURATION: 

ACCURACY: 

PRESET PARAMETERS: 



DESCRIPTION: 



METHOD: 



Irregular Spherical Bess el Functions (Spherical 

Neumann Functions) (DOI or SADOl) 

Closed, standard entry, with one program parameter 

29 

0, 1, 2, 3, h 

19.6 + 3.7N milliseconds 

7 decimal digits 

S3: during readin of this routine, location 3 must 

contain the address where the table of functions 

n (x), el (x), ...*;...- n^x) is to be stored. 
S^5 address of the first word of a closed sine-cosine 

subroutine which computes l/2 sin x (usually 

routine T5-157 will be used). 
To compute the Spherical Neumann Functions n„(x), 
n n (x), ... n^(x), enter this routine with the scaled 
argument 2 ^x in the accumulator and with the orders 

1 < N < 15 " 
< x < 32 (see below). 





— NF 


p 


50 pF 




26 qF 


p+1 


-- - 



where q is the location of this routine. When control 
is returned to the right hand side of p+1, n k ( x ) • 2 
will be found in location kS3. The functions n Q and n 
are always computed. 
The functions are computed using the formulas 

/ \ COS X 

(x) = - 



n, 







x 



cos x sin x 



, and 



n-^x) = £- 

x 

^+1 = x~~ "k " Vl f0r k > 1 ' 



- 2 - 

Because of the scaling used, all function: values 
must be less than kOS6 in absolute magnitude. 
Since n, (x) -> - co as x -» 0, x must be greater 
than the value given below in order to compute the 
functions up to a given N: 



If N is 


Then x must be greater than 


1 


0.01563 


2 


0.091 


3 


0.25 


k 


0.^9 


■5 


0.79 


6 


1.16 


7 


1.57 


8 


2.03 


9 


2*52 


10 


3.0^ 


11 


3.60 


12 


4,17 


13 


Ik 77 


Ik 


5.38 


15 


6.01 
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nj 



LOCATION 


ORDER NOTES PAGE 1 


V 16 




00K(Vl6) 











40 3F 
K5 IF 


by lV, 15* 


Store 2~ 5 x at 3 




1 


42 18L 
10 19F 




Plant link 
Nx2" 58 




2 


42 26L 

19 5F 




Plant 2N x 2~ 59 




3 


50 3F 










74 28L 








4 


00 5F 
40 4F 




(x/jt)(mod 2), rounded 




5 


LJ 4F 
50 5L 




jump to sin-cos subroutine 




6 


26 S4 




for 1/2 cos x 




7 


io i6f 
66 3F 
SI S3 


by 13' 


-17 
2 ' cos x ^-12 cos x 

o"5 . ~ x 
2 x 




8 


40 S3 

41 27L 




Store n Q (x) 

clear k count to zero 




9 


L5 4f 
50 9L 




jump to sin-cos subroutine 




10 


26 S4 
10 11F 




for l/2 sinx 




11 


LO S3 
10 5F 




-17r cos x , . . -> 
2 ' [ — — — H-rsmn xj 




12 


66 3F 


• -2. o~5 v 
! 1? 






SI F 








13 


40 1S3 
L5 7L 




store n, (x) 
set addresses 




Ik 


42 22L 
L4 l 




plant il address 




15 


42 20L 
L4 L 




plant n, address 




.16 


42 24L 
F5 27L 




plant n^ address 

step k count 





LOCATION 


ORDER 




NOTES PAGE 2 V 16 


17 


^0 27L 
Fk 27L 




■ 


18 


LO 26L 




2k+l-2N 




32 [ ]F 


by 1 


exit via link 


19 


IA 26L 








50 27L 




clear Q 


20 


10 5F 




(2k+l) 2" 5 in Q 


_ 


75 [ ]F 


by 15 


2 -17 (2k+l) n 


21 


40 F 




save MSP 




S5 F 




LSP to A 


22 


50 3F 








70 C ]F 


loylk 


-^\.i ] 2 " 17 


23 


L4 f 




(2fal) v --Tn tii 




66 3F 




X 


2k 


S5 F 








^0 [ ]F 


by 16 


Store n^ 


25 


F5 22L 




Advance addresses 




26 llj-L 




and loop 


26 


00 F 








00 [ ]F 


by 2 


2N x 2 -59 


27 


00 F 








00 F 




k x 2" 39 


28 


28 30^6F 








OJ 2964F 




l/it = .31830988618^ 



